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Abstract. We present astrometry and BVI photometry, down to V — 22, of the 
very young open cluster NGC6530, obtained from observations taken with the 
Wide Field Imager camera at the MPG/ESO 2.2 m Telescope. Both the V vs. 
+3 ■ B — V and the V vs. V — I color-magnitude diagrams (CMD) show the upper 

main sequence dominated by very bright cluster stars, while, due to the high ob- 
scuration of the giant molecular cloud surrounding the cluster, the blue envelopes 
of the diagrams at V > 14 are limited to the main sequence stars at the distance 
of NGC 6530. This particular structure of the NGC 6530 CMD allows us to con- 
clude that its distance is about d ~ 1250 pc, significantly lower than the previous 
determination of d=1800pc. 

We have positionally matched our optical catalog with the list of X-ray sources 
found in a Chandra- ACIS observation, finding a total of 828 common stars, 90% 
of which are pre- main sequence stars in NGC 6530. Using evolutionary tracks of 



Sicss et al 



2000), mass and age values are inferred for these stars. The median 
age of the cluster is about 2.3 Myr; in the mass range (0.6-4.0) Mq, the Initial 
Mass Function (IMF) shows a power law index x — 1.22 ± 0.17, consistent with 
both the Salpeter index (1.35), and with the index derived for other young clusters 
; towards smaller masses the IMF shows a peak and then it starts to decrease. 



Key words. Open clusters - individual: NGC 6530 - photometry - astrometry - 
membership - pre-main sequence stars - Luminosity and Initial Mass Function 
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1. Introduction 

Star formation regions and very young open clusters are crucial systems for understanding 
the star formation process because they allow us to derive Initial Mass Functions that are 
not affected by stellar and/or dynamical evolution effects. NGC 6530 (/ = 6.14, b = —1.38) 
is an example of a very young open clust er located in fron t of the M8 giant molecular 



cloud, also referred to as Lagoon Nebula l|Lada et al 



1976) . The brightest part of M8 



is the Hourglass Nebula, illuminated by the extremely young O star Herschel 36 (07 V). 
Other stars of spectral type O, 9 Sgr (04 V) and the binary system HD 165052 (06.5 
V + 06.5 V) excite the Lagoon Nebula. The open cluster NGC 6530 is located about 
9'.8 eastwards of the Hourglass Nebula (corresponding to a projected distance from the 
Nebula of 3.6 pc, using the cluster distance d = 1250 pc estimated in this work) and can 
be recognized by its brightest stars. 
Since the work of 



Walker (1957), several investigations have been devoted to study 
this cluster and to estimate its parameters. Using mainly photoelectric and photographic 
observations, many attempts have been done to estimate the distance of NGC 6530 from 
the CMDs, which show a normal cluster main sequence down to about AO stars; fainter 
stars lie above the main sequence, indicating th at the cluster i s so young that its low mass 



1957). The lack of the complete 



members are still gravitationally contracting (Walker, 
cluster main sequence made it hard to obtain reliable estimates of the cluster distance. 
NGC 6530 distance has been estimated by several authors to be in the range (1300- 
2000) pc (cfr. Tabled, by taking advantage of spectroscopic observations, to determine 
spectral types, and proper motions, to select cluster members. 

On the contrary, UBV observations have allowed several authors to study the red- 
dening law and to derive cluster redd ening values very similar to the average value 
E(B - V) = 0.35 recently derived by 
ground reddening value, E(B — V) f g 



Sung et al 



(2000), who have assumed a fore - 

(199oJ). 



0.17, as estimated by 



McCall et al 



Finally, the age of the cluster was estimated to be in t he range (1.5-2.0 ) Myr 



eral authors (van Altena fc Jones 



Damiani et al 



1972: 



Saear fc J oshi. 



1978: 



Sung et al 



by sev- 



2000), while 



2004) give a median age of 0.8 Myr. Due to the young age of this cluster, 



low mass stars are expected to be still in pre-main sequence phase, as already found in 



Sung et al 



I 2000J) by using photometric data down to the limiting magnitude V ~ 17. 



119 X-ray point s ources in the Lagoon Nebula region have been recently detected by 



Rauw et al 



I 2002) in a 20ksec XMM~Newton observation; they found that most of the 
X-ray sources are associated with pre-main sequence stars of low and intermediate mass. 
However, a larger list of point sources in the same region, with a much better spatial 
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Distance (pc) (m — M)o Reference 



1380-2000 


10.7-11.5 


Walker (1957) 


1300 


10.6 


The (19601 


1400 


10.7 


Walker (1961) 


1580 


11.0 


Hiltner et al. (1965) 


1780 


11.25 


van Altena & Jones (1972) 


1380 


10.7 


Kilambi (1977) 


1820 


11.3 


Sasrar & Joshi (1978) 


1900 


11.4 


Chini & Neckel (1981) 


1860 


11.35 


McCall et al. (1990) 


1600-2000 


11.0-11.5 


van den Ancker et al. (1997) 


1400 


10.74 


Loktin et al. ( 1997) 


1800 


11.25 


Sune et al. (2000) 


560-711 


8.75-9.26 


Loktin & Beshenov (2001) 



resolution, has been recently obtained by 



D amiani et al 



(2004), using Chandra ACIS-I 



X-ray data. Damiani et al. show that their source sample is made for at least 90% by 
cluste r members and tha t out of the 884 X-ray sources found, only 220 have a counterpart 



the 



Sung et al 



1 2000) optical catalog and most of this latter are pre-main sequence 
members of NGC 6530. As suggested by Damiani et al., the remaining X-ray sources are 
likely stars with magnitude fainter than V ~ 17. The lack of published photometric data 
at magnitudes fainter than V ~ 17 has motivated the analysis of deep optical images to 
study the low mass stellar population in the NGC 6530 field and, in particular, the low- 
mass pre-main sequence stars of the cluster by means of cross-correlation with available 
X-ray data. 

In this paper, we first present the observations and the data reduction procedure 
(Sect. 2) and the CMDs obtained using the optical catalog (Sect. 3). Therefore, we present 
the cross-correlation of optical and X-ray data, from which membership is determined 
(Sect. 4), and the cross-correlation of optical and 2MASS IR data (Sect. 5); mass and 
age determination of cluster members and the analysis of the spatial distribution are 
discussed in Sect. 6. Finally, the Luminosity and the Initial Mass Functions (Sect. 7) and 
our conclusions (Sect. 8) are presented. 



2. The Observational Data 



The data used in this work come from the combination of optical BVI images taken with 
the Wide Field Imager (WFI) camera at the 2.2 m Telescope of the European Sout hern 
Observatory (ESO), a 60 ks Chandra ACIS X-ray observation IjDamiani et all 120041) and 
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Table 2 Log-book of the optical observations. 



Target 


RA (J2000) 


Dec (J2000) 


Night 


Filter 


Exp. Time 


seeing 


Airmass 


EIS name 


(h m s) 


(d m s) 






[s] 


FWHM["] 




NGC 6530 


18 04 48.0 


-24 19 41.0 


27-28 Jul 2000 


B 


1 x 30 + 2 x 240 


1.36-1.76 


1.268-1.236 


(OC31) 


11 


)) 




V 


1 x 30 + 2 x 240 


1.19-1.36 


1.216-1.194 




)) 


)) 




I 


1 x 30 + 2 x 240 


0.83-1.05 


1.175-1.154 



public near infrared data f rom the All-Sky C atalog of Point Sources of the Two Micron 



All Sky Survey (2MASS) 1 ijCutri et all 120031 available on the WEB 2 



2.1. Optical Observations and Data Reduction 

The optical observations, consisting of 9 BVI images of NGC 6530, were taken using 
the WFI camera mounted at the Cassegrain focus of the ESO 2.2 m Telescope at La 
Silla (Chile). This instrument consists of a 4 x 2 mosaic of CCDs of 2048x4096 square 
pixels with a scale of 0.238 arcsec/pixel; each chip is 8' 12 x 16'25, while the full field 
of view (FOV) is 34 x 33 square arcmin. The observations, retrieved from the ESO/ST- 
ECF Science Archive Facility, ar e part of the ESO Ima ging Survey (EIS) in the context 



of the PRE-FLAMES program ( M omanv et al 



2001). Details of the observations are 



given in Table [2 The frames were taken under slightly variable seeing conditions with 
an improvement during the night of the Full Width Half Maximum (FWHM) of the 
point-spread function (PSF) from 1.76 to 0.83 arcsec, as measured on the frames. 

Fig. ^ shows the image obtained from the combination of the two deep dithered I 
band images of the region around NGC 6530, corrected for the instrumental signatures 
as described below. The FWHM of the PSF in this frame is about 0.83 arcsec. The 
very young cluster NGC 6530 can be recognized by its brightest stars located about 9'8 
eastwards of the O star Herschel 36, embedded in the Hourglass Nebula. The image shows 
regions of high and low stellar surface density according to the irregular pattern of the 
molecular cloud absorption. 

The first stage of the data reduction process has been the instrumental calibration of 
the images through the mscred package, a mosaic specific task implemented as an IRAF 3 
package for the NOACfi Mosaic Data Handling System. First the instrumental electronic 



1 a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics and Space 

Administration and the National Science Foundation 

2 http://irsa.ipac.caltech.edu/ 

3 Image Reduction and Analysis Facility 

4 National Optical Astronomy Observatory 
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Fig. 1 Image obtained from the combination of the two deep dithered / band images 
of the region around NGC 6530 observed with the WFI. White squares are unexposed 
regions due to the deterctor gaps. 



bias of the images was subtracted using the overscan region and then the images were 
trimmed to remove the overscan region. Flat fielding for each filter was performed using 
a set of sky flat fields scaled to the median value of all chips combined into a master flat 
using the mscred f latcombine task. 

A special treatment was required for the images in the / band, affected by strong 
effects of fringing. To remove this instrumental artifact, the fringing pattern, provided 
by the MPG/ESO 2.2 m Telescope team, 5 was subtracted scaling it to each exposure. 

Detection of sources recorded in the analyzed digital frame s was ob t ained using the 
DAOPHOT II/ALLSTAR photometric routines described in Istetsonl jl987jl . A more 



accurate PSF fitting photometry was performed by submitting all images to the routine 



5 available at http://www.ls.eso.org/lasilla/Telescopes/2p2T/E2p2M/WFI 
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ALLFRAME (Stetson 



1994f) that simultaneously uses the geometric and photometric 
information from all images to derive a self-consistent set of positions and magnitudes 
for all detected sources. 

In the present case, all images are characterized by a very strong sky background 
gradient and thus particular attention was devoted to estimate the minimum sky level 
(DAOPHOT parameter LOW GOOD DATUM=50 in unit of standard deviations) in order to 
include legitimate sky pixels where the background is faint. 

A non standard procedure, specifically suited for this case, was applied to the chip 
containing the star Herschel 36 because a significant number of spurious detections was 
found due to the emission from the nearby Hourglass Nebula. In order to recognize true 
star detections, the list of the stars released by ALLSTAR was filtered according to the 
sharp parameter. After this selection was applied, the simultaneous PSF fitting pho- 
tometry was performed with the routine ALLFRAME which releases a star-subtracted 
image from each frame. The star-subtracted images were used to obtain a median image, 
that was then smoothed with a Gaussian smoothing a larger than the FWHM of the 
PSF but smaller than the angular scale of the structures in the nebula. Finally a median 
image without nebula was obtained by subtracting the smoothed star-subtracted me- 
dian image from the median image obtained by the original images. A more reliable star 
list, determined using the latter median nebula-subtracted image, was given as input to 
ALLFRAME to obtain the final PSF fitting photometry. 

In order to convert the profile-fitting photometry to the standard photometric system, 
a magnitude zero point w as calculated f or each chip. This was done using the growth 



curve method described in 



Stetson! 1)19901) . First, we have selected the 5268 stars used to 



define the PSF model; next, all other objects were removed from the fram es and aperture 



photometry was carried out at different radii. The DAOGROW code l|Stetsonl . 11990 ) 
was used to derive growth curves and COLLECT was used to calculate the "aperture 
correction" coefficient for each chip, from the difference between PSF-fitting magnitudes 
and aperture photometry magnitudes of the selected stars. 

The procedure described above was also applied to a set of images of the 



Landolt 



( 1992) standard fields SA92 and SA107 obtained with WFI/2.2 m during the same night. 



Using the v, b and i instrumental magnitudes and the V, B and / magnitudes of the 
Johnson-Kron-Cousins photometric system of the standard stars falling in these fields, 
the transformation coefficients to the standard system were performed using the following 
equations: 



v =V + A + A 1 xX + A 2 x{ V-I), 
b=B + B + B 1 xX + B 2 x (B-V), 
i = I + Co + Ci x X + C 2 x ( V-I). 



(1) 
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where X is the airmass and Aq, Bq and Co are the magnitude zero points, A%, B\ and 
C\ are the extinction coefficients, while A2, B2 and C2 are the color terms. The list 
of standard stars used to determine the coefficients contains the standard stars defined 



the 



Landolt 



Landolt fields and defined by 



I i992) catalog and a number of secondary standard stars present in the 
Stetson (2000). Of these stars, we have between 25 and 
161 measurements for each chip, except for chip # 53 (the north-western chip in Fig. 
^ because neither Landolt's catalogue standards nor Stetson's catalogue standards were 
falling in this chip. Since the observations of the standard stars do not cover a wide 
airmass range, we were not able to derive the extinction coefficients. The extinction 
coefficients typical of La Silla (Stetson, private communication) were adopted, and the 
zero point and the color term were calculated for all chips except for chip # 53 . Stars 
with magnitude residual values above the 3 a level were not considered for the coefficient 
fitting. Considering the homogeneous characteristic of the instrument, we choose to com- 
pute new zero points fixing the color term to the mean value of the color terms of the 
chips # 50, #51, #54, #55, #56 for which the number of measured standard stars was 
larger than 40. Differences smaller than 0.1 mag among the zero points of different chips 
were found. Based on this latter result, we used their average and the above determined 
average color term to calibrate the stars in chip # 53, where no standard stars are found. 
The resulting coefficients are given in Table 



Table 3 Coefficients of the transformation to the standard system for each filter and 
averaged over all chips. 

Filter Av. Zero Point Extinction Av. Color Term 

V 0.905 ±0.021 0.14 0.066 ±0.046 

B 0.528 ±0.030 0.25 -0.245 ±0.019 

I 1.904 ±0.009 0.09 -0.133 ±0.034 



Using these coefficients and the inverted form of equations QJ, the B, V, /magnitudes 
of all objects detected in at least two filters were calculated. In order to test the self- 
consistency of the photometry, we compared the photometry of the stars common to 
contiguous chips. Common stars between chips were found because the deep mosaic 
images where dithered by about 30 arcsec in RA and about 60 arcsec in Dec from one 
another, in order to cover the inter-chip gaps. By matching contiguous chips we found 
about 10-30 common stars per chip pair for which we compared the photometry. We 
found mean offsets in the photometric zero points < 4% in V and /, and < 7% in B; this 
supports the choice to use the average zero point of the other chips to calibrate the data 
in chip #53. 
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In order to consider only stellar-like ob jec ts, a data selec tion was applied by including 



only objects with sharp between -0.5 and 0.5 (Stetson, 



1987j) i.e. objects with a brightness 



distribution consistent with a point source. After this selection, our catalog includes 
53 581 objects. 



2.1.1. Astrometry 

Celestial coordinates of the detected optical stars were obtained using the Guide Star 
Catalogue Version 2.2.01 (CSC 2.2, STScI, 2001) as reference catalogue. The first step 
was to match the list of celestial coordinates of stars retrieved from the GSC 2.2 catalog, 
with the list of pixel coordinates obtained for each chip as described in the previous 
section. The transformation between the two systems was obtained applying the appro- 
priate projection to the celestial coordinates and a linear transformation to the pixel 
coordinates. The initial estimate for the linear transformation is determined using three 
stars whose celestial and pixel coordinates are both known. To find the best astromctric 
solution, we matched the two catalogs using a conservative matching tolerance of 0"3. 
With this matching tolerance 2671 common stars were found (IRAF task ccxymatch) 
over a total of 33 554 sources (appearing in at least 50% of our frames) used to match the 
GSC 2.2 catalog. The matched lists of pixel and celestial coordinates of each chip were 
used to compute the plate solution. The sky projection was obtained using a combina- 
tion of the tangent plate projection and polynomials (IRAF task ccmap). Finally, using 
the derived astrometric solution, we computed celestial coordinates of our photometric 
catalog stars. 

In order to estimate our astrometric accuracy we matched stars of our catalog with 
V < 20, with the GSC 2.2 catalog, used as reference to find the astrometric solution 
and we considered the offset distribution within a relatively large value (5"). From this 
distribution we subtracted the expected distribution of spurious identifications, obtaining 
the distribution of true identifications only. The resulting rms is 0"4, that is therefore 
our final accuracy. 



2.1.2. Comparison with previous catalogs 

As an external check of our photometry we sel ected stars of our catalog with V < 17.5 



and compared them with the 



Sung et al 



2000) catalog limited to V < 17. By matching 
the two catalogs, we found a systematic offset in both right ascension and declination of 
— 1"05 and ~0"34, respectively. 

To estimate the appropriate matching radius for comparing the two catalogs, we first 
removed the coordinate offsets and then considered the offset distribution within 7". We 
found that the rms width of the distribution is 0"5 and that the contribution of spurious 
identifications becomes dominant at offsets larger than 1". 



Prisinzano et al.: NGC6530 



With the aim to include only reliable matches between our and the Sung et al. cat- 
alogs, we adopted a matching radius of 0"5, finding 576 stars with offsets -0"01±0"30 
and -0''02±0"22, in RA and Dec, respectively. The mean offsets and the rms of the pho- 
tometric residuals are 0.006 ± 0.060 mag, 0.007 ± 0.075 mag and 0.037 ± 0.055 mag in 
V, (B — V), and (V — I), respectively. These values are computed after applying a 3<r 
clipping to the photometric residuals. 



Final l y, we 



Kilamb 



compared our photometric catalog with the catalog of 



Walker (1957) and 



1 1977) consisting of a total of 319 stars. Of these stars, only the 150 stars with 
V > 9.5, falling in the WFI FOV and with celestial coordinates found and retrieved from 
the SIMBAD database, were considered. With a matching radius of 9" we found 131 
stars plus a number of mis-identifications that we discarded based on the ir photometric 
resi duals. By consid ering the stars having both B and V magnitudes in the 



and 



Kilamb 



Walker ( 1957) 



19771) catalog, we found = 0.002 ±0.073 and A(B-V) = 0.006 ±0.081, 



after applying a 5<r clipping to the photometric residuals. The remaining 19 stars with 
no counterpart in our catalog have V ~ 10, hence they are saturated in the WFI frames. 



All previous comparisons show that our magnitudes and colors are in excellent agree- 
ment with those measured by other authors. Table 01 shows the results of such compar- 
isons; col.s 1-9 list the sequenti al number, celesti al coordinates and photometry of our 



catalog; the star number used by 



Sung et al 



1 2000, denoted as IDscr) and the astromet 



ric and photometric residuals obtained from the compariso n with th e 
catalo g are given in col.s 10-15; the star number used by 



Sung et a,iy.2000) 



Walker (1957) and Kilamb 



1 19771 denoted as IDwk ) and th e astro metr ic and photome tric residuals obtained from 



the comparison with the 
16-20. 



Walker (1957) and 



Kilamb 



(197JJ) catalogs are given in col.s 



2.2. X-ray and near-infrared data 



To better underst a nd pa rt of our optical data we used a list of X-ray sources, presented 



D am iani et al 



I 2004) . These data were obtained from a 60 ks observation carried 
out on June 18-19, 2001 with the Chandra ACIS-I CCD detector having a very narrow 
PSF (FWHM-0.5" on-axis). The observation consists of a FOV of 17'xl7' pointed 
toward RA=18' l 04 m 2Ti38, Dec=-24°21 / 05'?8, i.e. centered on t he NGC 6530 cluster. 
point-like X-ray sources were detected as described in detail by 



Damiani et al 



(2004) 



In addition, to complement the optical data, we used JHK$ infrared (IR) photometry 
taken fro m the All- Sky Point Source public catalogue of the Two Micron All Sky Survey 
(2MASS) icutri et aljbflj). 



o 



Table 4. Cross-i dentif i cations of this c atalo g with the previou s works of lWalker 



(1957); 



Kilambil 1 1977) and 



Sung et al 



(2000) 



RA(2000) Dcc(2000) 
(h m s) (dm s) 



A RA a A Dec" Ay* A(S - V) a A(V - J) a ID WK A RA b A Dcc b 
(arcscc) (arcsec) (arcscc) (arcscc) 



AV A(S - V)" 



WPI 


5173 


18 


04 


35 


45 


-24 


32 


28.82 


14 


74!) 





002 


1.3 


889 





003 


12 


800 





004 


732 


0.29982 


-0 


36392 


-0.061 


0.016 


0.022 


WFI 


5249 


18 


04 


45 


36 


-24 


32 


26.44 


17 


692 





005 


10 


094 





003 


15 


477 





006 


836 


0.39976 


-0 


69351 


-0.028 


-0.040 


0.006 


WFI 


5377 


18 


04 


33 


15 


-24 


.32 


22.37 


18 


709 





004 


10 


001 





003 


14 


008 





003 


712 


0.39976 


-0 


39825 


-0.036 


0.177 


0.072 


WFI 


5387 


18 


04 


13 


26 


-24 


32 


22.09 


15 


973 





000 


15 


115 





006 


14 


091 





003 


451 


0.49970 


-0 


34332 


-0.108 


-0.007 


0.009 


WFI 


5449 


18 


04 


31 


27 


-24 


32 


20.17 


13 


460 





008 


13 


049 





009 


12 


605 





012 


087 


0.69958 


-0 


32959 


-0.038 





009 


0.030 


WFI 


5472 


18 


03 


47 


86 


-24 


32 


19.58 


17 


853 





00.3 


10 


789 





003 


15 


587 





006 


207 


0.29982 


-0 


52872 


-0.080 





028 


0.024 


WFI 


5483 


18 


04 


37 


62 


-24 


32 


19.21 


15 


493 





02.3 


14 


305 





012 


12 


811 





031 


750 


-0.19988 


-0 


26779 


0.001 





181 


-0.021 


WFI 


5493 


18 


03 


51 


58 


-24 


.32 


18.72 


15 


619 





00.3 


14 


779 





001 


1.3 


778 





004 


247 


-0.69958 


-0 


46005 


-0.061 





0.35 


0.024 


WFI 


5556 


18 


03 


48 


69 


-24 


32 


16.89 


15 


556 





002 


14 


596 





003 


1.3 


.327 





005 


218 


-0.49971 


-0 


53558 


-0.107 





124 


0.030 


WFI 


5660 


18 


03 


48 


32 


-24 


32 


13.42 


14 


556 





005 


1.3 


645 





003 


12 


635 





010 


212 


-0.19988 


-0 


44632 


-0.069 





033 


0.017 


WFI 


5828 


18 


03 


43 


17 


-24 


32 


7.96 


18 


018 





00.3 


10 


018 





005 


1.3 


081 





006 


168 


-0.69960 


-0 


50125 


-0.031 





129 


0.102 


WFI 


5857 


18 


03 


36 


99 


-24 


32 


7.15 


17 


187 





00.3 


14 


814 





004 


11 


801 





010 


117 


-0.49972 


-0 


47379 


-0.005 





14.3 


0.129 


WFI 


5875 


18 


03 


50 


34 


-24 


.32 


6.68 


19 


092 





007 


10 


951 





005 


14 


451 





003 


234 


0.00000 


-0 


51498 


0.021 





091 


0.114 


WFI 


5965 


18 


03 


45 


52 


-24 32 3.92 


18 


660 





004 


10 


005 





003 


14 


090 





005 


187 


-0.09994 


-0 


45319 


-0.049 





128 


0.083 


WFI 


5979 


18 


03 


30 


39 


-24 32 3.45 


17 


355 





011 


10 


313 





007 


15 


024 





010 


83 


-0.59966 


-0 


70038 


-0.141 





045 


0.031 


WFI 


5995 


18 


03 


51 


04 


-24 


32 3.01 


17 


879 





005 


15 


410 





002 


12 


450 





013 


241 


0.19989 


-0 


44632 


0.028 





122 


0.102 


WFI 


6004 


18 


03 


511 


02 


-24 


.32 


2.75 


16 


633 





002 


15 


002 





002 


14 


400 





003 


303 


0.00000 


-0 


40512 


-0.034 


-0.027 


0.037 


WFI 


6057 


18 


04 


2!) 


68 


-24 32 0.99 


19 


260 





000 


10 


04.3 





004 


1.3 


365 





004 


075 


0.79956 


-0 


43945 


0.071 


0.072 


0.140 


WFI 


6149 


18 


04 


20 


02 


-24 31 


57.73 


17 


911 





002 


10 


289 





006 


14 


202 





002 


541 


-0.29984 


-0 


45319 


-0.005 


0.004 


0.058 


WFI 


6278 


18 


03 


43 


40 


-24 31 


52.64 


18 


538 





00.3 


10 


515 





004 


14 


071 





005 


172 


-0.19989 


-0 


47379 


-0.049 


0.064 


0.077 



3.09780 -2.10980 -0.031 



N 

ts 
o 

O 
O 
o 
^> 
w 
o 



a SCB comparison 

b Walkcr-Kilambi comparison 

Note. — The complete table is available in electronic format at the Web page http://cdswcb.u-strasbg.fr/ 
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3. The Color-Magnitude Diagrams 

3.1. Optical data 

The fi nal optical photom etric catalog has been obtained by adding to our data, the stars 



of the 



Sung et al 



I 2000) catalog that were not detected in the WFI images. These are 
very bright stars (V < 10 mag), saturated in our images, or stars falling out of the 
WFI FOV or on the edges of t he chips, fo r a tot al of 123 objects. To select these stars, 



we considered the stars of the 



Sung et al 



I 2000) catalog whose astrometric offset with 



respect to our catalog is larger than 1", the radius for which the number of matched stars 
drops to zero. 

CMDs of all stars of the complete optical catalog were constructed after selecting 
the stars with errors in the V and / magnitudes both smaller than 0.2 mag; stars with 
larger (statistical) errors cannot be placed accurately on the CMD so their inclusion 
is meaningless. We applied this selection only to the V and I magnitudes because in 
the following analysis we mainly will use the V vs. V-I CMD. Nevertheless, the V vs. 
B-V diagram has been obtained using stars of the cleaned catalog with errors in B 
smaller than 0.2 mag. Using this selection criterion on our data we have found that our 
detection limit is V ~ 23. Completeness of our data was determined via artificial star 
tests. A total of about 1650 artificial stars were added on each chip and photometry 
of the artificial frames was performed using the same reduction procedure applied to 
the original frames. In order to test the completeness of the data used in the following 
analysis, the results were obtained by considering the number of retrieved artificial stars 
with sharp between -0.5 and 0.5 and errors in the V and / magnitudes smaller than 0.2 
mag. We have found that our data are 100% complete above V ~ 20 and / ~ 18.5, while 
they are more than 50% complete above V ~ 22 and I ~ 19.5. 

The V vs. B-V and the V vs. V-I CMDs, obtained using the cleaned catalogs, are 
shown in Fig. where horizontal bars indicate the median errors in color, while vertical 
bars (barely visible), indicate the median errors in magnitude for one magnitude bins. 
The number of stars and the error bars in the V vs. B-V diagram are smaller than those 
in the V vs. V-I diagram because in the V vs. B-V the selection has been applied to 
the B, V and I magnitudes while in the V vs. V-I diagram it has been applied only to 
the V and I magnitudes. 

The hor izontal dashed line in each diagram, indicating the magnitude completeness 



limit of the 



Sung et al 



( 2000) catalog, gives an idea of the amount of new photometric 
information yielded by this survey. 

Both diagrams clearly show the upper main sequence (B-V< 0.5 and V-I< 0.6) 
dominated by the very bright cluster stars. These are the more massive stars which, 
having a very short pre-main sequence lifetime, have already reached the main sequence 
showing a very small age spread. 




Fig. 2 The V vs. B - V and the V vs. V - I CMDs of selected stars (20 253 and 37 553, 
respectively) in the WFI FOV. Horizontal bars indicate the median errors in color, while 
vertical bars (barely visible) indicate the median errors in magnitude for bins of one 
magnitude. The horiz ontal dashe d line in each diagram indicates the magnitude com- 

2000) catalog. The solid line is the Zero Age Main 



pleteness limit of the_ 
Sequence (ZAMS) of 



Sun g et a 



Siess et al 



I 2000) assuming the cluster distance d — 1250 pc and 



E(B — V) = 0.20; the thick long dashed line is the same ZAMS at th e same distance bu t 
with the average cluster reddening E {B — V) = 0.35 estim ated bv ISung et alJ (2000); 
finally the dashed line is the ZAMS of Sung fc BesselJ 1 1999) at the distance of 1800 pc 
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As already noted by 



Sung et al 



2000), a sequence of stars ap pears in the 



Sung et al 



range 



(2000) 



0-7 <V — 7 < 1.5 and 13 < V < 18. These stars were interpreted by 
as foreground stars less reddened than the remaining stars in the field. 

Cluster stars fainter than V ~ 14 are expected to be pre-main sequence stars which 
populate a large region in the CMDs, highly contaminated by field stars. 

One of the most important features of these diagrams is the well defined blue envelope; 
it is due to the presence of the giant molecular cloud, which prevents to see field stars 
(mostly main-sequence) more distant than the cloud. Therefore, the well defined blue 
envelope of the CMDs is populated by main sequence field stars at the distance of the 
cloud. Background field stars, i.e. stars more distant than the cloud, are highly obscured 
by the cloud and therefore they would be visible at magnitudes and colors much fainter 
and much redder than their intrinsic values. 

Finally, we note the presence of an unexpected sequence of red stars, well separated 
from the bulk of the other objects. The analysis of the spatial distribution of these 
stars indicates that they are uncorrelated with the region dominated by the cluster stars, 
identified by a clustering of selected members, as will be shown in Sect. El On the contrary, 
they cover almost uniformly the less obscured regions of the molecular cloud and are likely 
associated with clump red giants at larger distances. Since we are here mainly interested 
to study the stellar population of the very young cluster NGC 6530, we defer a further 
investigation of these stars to a future work. 



3.1.1. Cluster distance 

As already mentioned in the previous section, at NGC 6530 age and distance, only stars 
brighter than V ~ 14 define the main sequence that in this magnitude range is al- 
most vertical preventing the estimate of the cluster distance using the usual photometric 
method. 

Fainter cluster stars do not define any sequence, but, on the contrary, they populate 
a large region of the CMD which cannot be used to constrain the cluster distance. 

Instead, we derived a reliable estimate of the cluster distance by taking advantage 
of the particular feature of its CMDs: the well defined blue envelope. It is populated 
by main sequence field stars at the cloud distance, which is very similar to the cluster 
distance if one does the very reasonable assumption that the cluster lies just before 
the cloud. In fact, the blue envelope fixes a magnitude limit within which we see stars 
(either belonging to the cluster or not), with a reddening equal or less than the cluster 
reddening. In particular, stars along the blue envelope brighter than V ~ 14 are cluster 
stars that have reached the main sequence, while fainter stars along the blue envelope 
are foreground stars up to the cloud distance. 
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Siess et al 



( 2000) on CMDs (solid 



We superposed the theoretical ZAMS computed by 
line in Fig.PJJ) and we found that the best match between this curve and the blue envelope 
of the star distribution is found if a distance of d ~ 1250 pc is adopted, corresponding to 
a distance modulus (V — My)o — 10.48, and a reddening E(B — V) = 0.20 (corresponding 
to Ay = 0.6 and E(V - I) = 0.26 using the reddening law E(V - I) = 1.3 x E(B - V) 



by 



Munari & Carrarol p996l^ . We note that the reddening E(B - V) = 0.20 is the value 



derived for foreground field stars, which suffer from a reddening significantly smaller than 
cluster stars. In fact, for V < 14, the theoretical model is bluer than the bright stars 
because these are main sequence cluster stars at the same distance (d ~ 1250 pc) but 
with reddening E(B - V) = 0.35 (correspon ding to A v = 1.1 and E(V - I) = 0.46), 



that is the average cluster reddening given by 



Sung et al 



I 2000), as shown by the thick 



long dashed line in F ig. HI For compari son the ZAMS at the distance of 1800 pc, i.e. the 



Sung et al 



I 2000), is shown as dashed line] the extinction vectors 



distance adopted by 

Ay = 0.6 and Ay = 1-1, corresponding to foreground stars and cluster stars respectively, 
are also shown in Fig. [21 

The distance value d ~ 1250 pc derived for NGC6530 is very close to the lower 
limit of the ample range of values (1300-2000 pc) photometrically determined for bright 
stars in previous works; i t is, instead, significantly larger than the value obtained by 



Loktin & Beshenovl 1)20011) from the Hipparcos trigonometric parallaxes of 7 stars (see 
Table UJ. 



4. Optical Data and X-Ray Source Identifications 



As already discussed in the Introduction, X-ray observations are v ery useful tra c ers of 
membership in star f ormation regi o ns or very young clust e rs (e. g. 



Randich et al 



1995 



Sterzik et al 



1995 



Flaccomio et al 



Alcal a et al 



1995 



2000), such as NGC6530. 



Alternative membership assessment methods are not suitable for NGC 6530, since neither 
reliable proper motion measures nor radial velocities for magnitudes fainter than V ~ 14 
are yet available. 

To trace the pre-main sequence locus of NGC 6 530, the optical data derived in this 



Damiani et al 



(2004), although the 



work were matched with X-ray sources published by 
Chandra ACIS FOV is about 4 times smaller than the WFI FOV. However, the Chandra 
ACIS FO V is centered very close to the NGC 6530 nominal center RA=18 h 04 m 25 s , Dec=- 



24°22'00"(Damiani et al 



2004), thus allowing us to study most of the cluster stars. 



The total number of o ptical sources falling in the Chandra ACIS FOV is 8956, while 



the 



Damiani et al 



1 2004) catalog contains 884 X-ray s ources. In orde r to es timate the 

(2004) examined 



number of X-ray sources that are not cluster members, 



Damiani et al 



an observation of the Galactic Plane obtained with the same detector and used it as 
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" control field" . By analyzing the detected count-rate distributions of the two fields they 
conclude that at least 90% of the X-ray sources are very probable cluster members. 

To cross-correlate the X-ray and optical catalogs we have used a matching distance 
d < 4<tx (where ax is the X-ray position error, see Damiani et al. 2004). In doing this, 
we found and corrected a systematic shift between X-ray and WFI positions of 0''2 in 
RA, and -0"26 in Dec. For three X-ray sources with such a small X-ray error that 4ax 
was less than 1"5, we relaxed the identification condition to d <1"5. This resulted in a 
number of multiple identifications, among which four turned into single identifications 
by using a reduced distance d < 1.5ax- This leaves us with 721 single, 44 double and 3 
triple X-ray identifications in the optical catalog; in addition, 1 X-ray source corresponds 
to 4 optical identifications and 1 X-ray source corresponds to 6 optical identifications. 
These multiple identifications are near the edge of the ACIS FOV where the spatial 
resolution is much worse than in the center. The total number of X-ray sources with WFI 
counterpart (s) is therefore 770; of them only 15 X-ray identified stars come from the Sung 
et al. catalog and are not in the WFI catalog. The total number of optical sources with an 
X-ray counterpart is 828. The agreement between X-ray and WFI positions is excellent 
in most cases, with offsets below 1". The final list, comprising 770 X-ray sources, with 
the corresponding optical counterpart (s), is given in Table [5p, where col.s 1 and 2 are 
the celestial coordinates, col. 3 indicates the parent optical catalog, col. 4 is the optical 
identification number, col.s 5-10 are the B VI magnitudes and th eir uncertainties, col. 



Damiani et al 



1 2004) catalog and, finally, 



11 is the X-ray identification number of the 
col. 12 is the number of optical sources corresponding to the X-ray detections. 

Fig. shows the V vs. V-I CMD of stars falling in the Chandra ACIS FOV. Dots 
indicate all stars in the optical catalog falling in the ACIS FOV, while large f illed symbols 



Sung et al 



indicate optical stars having an X-ray counterpart; stars added from the 
I 2000) catalog are marked by squares. As expected, X-ray detected stars trace very well 
the CMD region occupied by the cluster stars. The cluster region identifie d by 220 X-ray 



source s, at magnitudes brighter than V 
I 2004), using the photometric catalog of 



17, was 



Sung et al 



already analyzed in 



Dam iani et, al 



12000). With the present survey it 



is now possible to identify the fainter cluster pre-main sequence region. In fact, most of 
the stars identified with an X-ray source are located within a well defined region, whose 
width may be due to an age spread and/or to binary stars. Only a small percentage 
of stars with an X-ray counterpart has photometry that is not consistent with that of 
the cluster. The field stars contamination will be quantified in Section [SI but at least 
qualitatively, it a ppears consistent with the percentage of contaminating X-ray sources 



(< 10%) given by 



Damiani et al 



available in electronic format at the Web page http://cdsweb.u-strasbg.fr/ 
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Table 5 Cross-identifications of the optical catalog with the X-ray source catalog of 
Damiani et al. (2004). 



RA(2000) 
(h m s) 


Dec(2000) 
(d m s) 


Cat 


ID 


B 




V 




I 




ID(X) 


Mult. 

# 


18 4 48.38 


-24 29 3.94 


WFI 


10889 


17.346 


0.008 


15.901 


0.008 


14.011 


0.005 


811 


1 


18 4 26.49 


-24 29 0.30 


WFI 


10947 


20.717 


0.109 


18.975 


0.024 


16.096 


0.031 


528 


2 


18 4 26.46 


-24 28 59.95 


WFI 


10956 


18.314 


0.019 


16.729 


0.011 


14.771 


0.007 


528 


2 


18 4 30.41 


-24 28 53.16 


WFI 


11091 


20.006 


0.010 


18.289 


0.018 


15.930 


0.004 


609 


1 


18 4 26.29 


-24 28 48.61 


WFI 


11173 


20.225 


0.013 


18.809 


0.011 


16.370 


0.004 


524 


1 


18 4 21.80 


-24 28 46.25 


WFI 


11216 


22.255 


0.065 


20.345 


0.024 


16.760 


0.007 


417 


1 


18 4 31.37 


-24 28 44.49 


WFI 


11246 


21.783 


0.046 


19.634 


0.034 


16.536 


0.004 


628 


1 


18 4 22.54 


-24 28 40.75 


WFI 


11308 


15.464 


0.003 


14.705 


0.004 


13.693 


0.005 


432 


1 


18 4 16.76 


-24 28 36.94 


WFI 


11376 


11.803 


0.038 


11.580 


0.003 


11.102 


0.005 


305 


1 


18 4 23.34 


-24 28 31.98 


WFI 


11471 


19.082 


0.018 


17.726 


0.004 


15.938 


0.004 


451 


1 



As it will be shown in Section El the cluster pre-main sequence region is comprised 
between the isochrones of 0.3 and 10 Myr; we used these isochrones to select possible 
cluster members from our optical catalog. 



5. Infrared data 

As discussed in the previous section, background stars are heavily obscured by the molec- 
ular cloud and therefore they appear at magnitudes much fainter than their intrinsic 
values. Nevertheless, optical diagrams presented in the previous section do not allow us 
to distinguish very faint objects with no or very low reddening, from intrinsically bright 
background objects appearing as very faint stars due to their high reddening. Since near- 
infrared bands are very sensitive to reddened stars, we constructed the near-infrared 
CMDs using the JHKs magnitudes from the 2MASS public catalog, with the aim to 
better characterize the stellar populations in the WFI FOV. The diagram H vs. H — K 
of stars of the 2MASS catalog falling in the WFI FOV is shown in Fig. [J] We note 
that only stars with magnitudes measured from point spread-function fitting or aperture 
photometry were selected from the 2MASS catalog. 

The large reddening effect can be quantified by inspecting the H vs. H — K diagram 
In fact , in this diagram most stars are intrinsically colorless (0 < H—K < 0.5, 



Alves et al 



1 1998)) and therefore the apparent H — K color of a star is, to first order, a qualitative 
measure of the amount of extinction toward this star. This means that, apart from 
relatively few stars having reddening of the order of the average cluster reddening E(B — 
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Fig. 3 V vs. V - I CMD of stars within the Chandra ACIS FOV. Dots indicate all stars 
in the optical catalog falling in the ACIS FOV, while large fille d symbols ind icate optical 



Sung et al 



stars having an X-ray counterpart; stars added from the 
marked by squares. Lines and extinction vectors are as in FigEl 



2000ft catalog are 



V) ~ 0.35, corresp onding t o E(H — K) ~ 0.07, the bulk of the stars have colors H — K > 



0.3 that, using the iMathisI l)l990|) reddening law, corresponds to Ay > 4.4. 

In this diagram we superposed three theoretical curves, that are the ZAMS, the 0.3 
and 10.0 Myr isochrones of Siess et al., 200 (see Sect . Et reddened using the cluster 
reddening value E{B — V^)=0.35 and the iMathia l)199Cf) reddening law. The location of 



these curves can be used to identify cluster stars and foreground field stars, affected by 
a reddening equal to or smaller than the average cluster reddening, from the bulk of the 
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Fig. 4 H vs. H - K CMD of all 2MASS point sources falling in the WFI FOV. The solid 
lines are the isochrones of 0.3 and 10.0 Myr, whereas the dashed line is the ZAMS at the 
cluster distance d — 1250 pc. 



high reddened stars (H — K > 0.3) that are either stars beyond the molecular cloud or 
cluster members showing IR excesses, i.e. Classical T Tauri stars. 



5.0.2. Optical and 2MASS detections 

In order to better characterize the stellar population of the optical catalog, we matched 
optical data with the 2MASS catalog, using a matching radius of 0"8. We found system- 
atic offsets in both right ascension and declination of -0"3 and 0''16, respectively; after 
correcting the optical catalog for this offset, a total of 15 282 common stars were found. 

Fig. El and Fig. H3 are the color-magnitude and the color-color diagrams of possible 
cluster stars detected both in the optical and in the near- infrared catalogs. They are the 
9613 stars (dots) that, in the (V — I) vs. V diagram, are found between the isochrones 
at 0.3 and 10.0 Myr, that is the estimated age spread of the cluster (see Sect. EJ. Those 
stars for wh ich an X-ray counterp art was also found are indicated as black bullets and, as 



asserted in 



Damiani ct al 



12004), more than 90% of them are cluster members. In these 
diagrams, we also plotted the isochrones of 0.3 and 10 Myr, indicated by the solid lines, 
and the redde ning vector. The do tted line in Fig. corresponds to the locus of classical 



T Tauri stars IjMever et al 



1997). 
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If we compare Fig.s0]and[S]we note that most of the reddened stars (e.g. H — K > 0.6) 
in the IR catalog are not seen in the optical bands, at least within the limiting magnitude 
of this survey. Most of the stars detected both in the optical and in the near-infrared bands 
have H — K < 0.6, b ut there are objects with (H — K) up to 1.5. Using the extinction 



law of 



Mathii 



1 1990), these values indicate that most of the stars in the optical catalog 
have extinction Ay < 8.8 mag, while the high reddening objects can have Ay up to 20 
mag. 

Fig. [3] clearly shows that, many cluster members (black bullets) are located near the 
young theoretical isochrone, but a significant number is also found redwards, with (H—K) 
up to 1.3. Therefore, whereas the cluster members (selected for their X-ray emission) are 
located in a well defined region of the optical CMD, they are much more "dispersed" in 
the H vs. H — K diagram. This is probably due to the fact that a large fraction of cluster 
members shows significant excesses in the IR colors. 

This can be seen also in the H — K vs. I — J diagram of Fig. H3 where many of the 
2MASS-optical candidate members are found in a strip around the theoretical curves, 
while at low reddening a clump of stars is found at H — K > 0.4 and (I — J) > 2.5. As 
already discussed, the latter objects are the highest reddened stars and have (V — I) > 3.5 
in the V vs. (V — I) diagram. Among the 2MASS-optical sources that are also X-ray 
detected, less than 10% are high reddening sources (7 — J > 2.5). Many of the X-ray 
detected cluster members are instead found along the theoretical isochrones, but there 
is also a significant number of stars with large H — K, but in a different direction with 
respect to the reddening vector. This supports the hypothesis th at these X-ray source s 



Damiani et al 



(2004). 



are Classical T Tauri stars with IR excesses, as also suggested by 

These properties are found also in the other combined optical-IR color-magnitude 
and color-color diagrams and therefore we conclude that optical and 2MASS combined 
data alone cannot be used to select clean sample of PMS cluster members. 



6. Masses, Ages and Spatial Distribution 

Using the 



Sung et al 



Dam iani et al 



1 2004) determined mass 



1 2000) photometric catalog, 
and age of X-ray detected stars brighter than V = 17, finding stars down to 0.5M©, with 
completeness only for masses larger than 1.5M Q . In addition, they found that the median 
age of the cluster stars in the central region is 0.82 Myr with a maximum spread of 4 Myr. 
These values were derived assuming a cluster distance larger than the value estimated in 
this work, which led those authors to overestimate luminosity and mass of the stars. 

In the present work we have recomputed mass and age of the probable cluster 
members, derived by interpolating the theoretical tracks and isochrones calculated by 



Siess et al 



1 2000j) to the positions of the stars in the V vs. V—I CMD for which an 



X-ray counterpart has been found. Of these stars more than 90% are cluster members, 
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Fig. 5 H vs. H — K CMD of the stars detected both in the optical and in the near- 
infrared catalogs. Dots are stars with age between 0.3 and 10 million years as selected in 
the (V — I) vs. V diagram while black bullets are those with an X-ray counterpart. The 
solid lines are the isochrones of 0.3 and 10.0 Myr. 



although it is unknown which percentage they are of the whole cluster population. The 
completeness of this sample is discussed in Section 

We used the models of 



Siess et al 



(2000) because it was shown tha t masses pre- 



dicted with these tracks are consistent with masses estimated dynamically IjSimon et al 



2000). In additi on the y are in good agreement with the pre- main sequence models 



of 



Baraffe et al 



11998), which are available only for ages older than lMyr. We are, 
however, aware of the large uncertainties of these models for age younger than 1 Myr , 
which are based on oversimplified initial conditions, as discussed in lBaraffe et alJ (2002). 
Nevertheless, they are the most complete set of theoretical tracks avail able in the litera- 
ture, already used in very recent works on star formation regions (e.g. 



Flaccomio et al 



2003). 



Since metallicity has never been estimated for NGC 6530, we assume for this cluster a 
solar metallicity and we consider the 



Siess et al 



(2000) models with Z = 0.02, Y = 0.277, 



X = 0.703 and no overshooting. To convert effective temperatures and luminosities of 
the adopted model to the e mpirical V vs. V-I plane, we used the conversion table of 



Kenvon & Hartmannl (J.1995) where both optical and infrared colors are available. 
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H-K 
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Fig. 6 / — J vs. H — K color-color diagram of the stars detected both in the optical and 
in the near-infrared catalogs. Symbols are as in Fig.[5j The solid lines are the isochrones 
of 0.3 and 10.0 Myr, whe reas the dotted line corresponds to the locus of classical T Tauri 
stars JlVlever et all ll£ 



Fig.0shows the V vs. V-I CMD of X-ray detected cluster members with overimposed 
the adopted theoretical curv es. We note t hat th e curves were reddened using the value 
E[B - V) 



Sung et al 



1 2000) and in agreement with our data (see 



0.35 given by 

Sect . l3~Tll. Ay = 3.1 x E(B - V) and the reddening law E(V - I) = 1.3 x E{B 



V) <|Munari fc Carraro 



1996). The absolute magnitudes of the theoretical models were 



transformed to apparent magnitudes assuming the distance d ~ 1250 pc, derived in Sect. 
13.11 Dotted lines are the evolutionary tracks for masses between 0.25 and 7 Mq, while 
solid lines are the isochrones for ages between 0.1 and 100 Myr, this latter being very 
similar to the ZAMS (dashed line). 

Due to the strong variation of the surface spatial density of this area, the cluster is 
expected to be affec ted by a non-negli gible differential reddening. Evidences of such ef- 



fect were found by 



Sung et al 



I 2000). This implies that the age and mass estimates 



are affected by unce rtainties due to differential reddening. However, as discussed in 



Dam iani et al 



I 2004), for most of the low mass cluster members, the reddening vec- 
tor direction is almost parallel to isochrones and thus the age estimates are less affected 
than mass estimates that can have an uncertainty up to 0.5 M@. Nevertheless, this un- 
certainty does not drastically affect the determination of the Initial Mass Function since 
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Fig. 7 V vs. V - I CMD of the stars within the Chandra ACIS FOV with an X-ray 
counterpar t . The dotted lines are the solar metallicity evolutionary tracks computed by 



Siess et al 



2000); the corresponding mass value in solar mass units is given at the redder 
extreme of each track. The solid lines are the corresponding isochrones for age between 
0.1 and 100 Myr, while the dashed line is the ZAMS at the distance of the cluster. We 
no te that the t heore tical curves were reddene d using the cluster redden ing value given 



by 



Sung et al 



1 2000j) and the reddening law of 



Munari fc Carrard (1996). 



it is smaller than the considered mass bins. Fig. [7| shows that most of the X-ray detected 
stars are included between the theoretical isochrones of 0.3 and 10 Myr. We expect that 
stars outside this age range are mainly contaminating objects or highly reddened cluster 
members. We note that the observed spread (1.5-2 mag) is larger than that expected 
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from binarity or photometric uncertainty since only data with errors smaller than 0.2 
mag have been considered. The age distribution for the stars in this age range is shown 
in Fig. El it is quite symmetric and concentrated, with a median age of 2.3 Myr. 
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Fig. 8 Age distribution of all optical/X-ray detected sources; the star age was estimated 
by comparing t he sta r positions in the V vs. V — I CMD, with the isochrones computed 



by 



Siess et al 



(2000). The median age value is 2.3 Myr. 



Due to the position of the cluster with respect to the Hourgl ass nebula, i t was 



posed for this cluster a sequenti al process of the star formation IjLada et all Il97fill . In 



particular, 



Damiani et al 



pro- 



1 2004|) suggested that the star formation has progressed from 
north to south. In order to find further evidences of this effect, we plotted in Fig. El the 
spatial distributions of the optical/X- ray detected sources i n four different age ranges. 



Our results support the conclusion of 



Dam iani et al 



since very young stars are 



present (almost) only in the cluster southern region, while the spatial distribution be- 
come more uniform for older ages. To test the validity of this conclusion, we have done a 
statistical analysis using the age distributions in 5 different spatial regio ns, namely the 
cluste r center included in a square region of side 7.2 arcmin, as defined in lDamiani et al. 
I 2004), and 4 outer regions (see Fig.^J. The cumulative age distributions in Fig. 1 101 show 
a clear age difference between the southern and central regions and the northern ones. 
Using the two sample Kolmogorov-Smirnov tests, we find that the cumulative distribu- 
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tion functions of each of the southern and central regions are different from the northern 
regions at a significance level greater than 99%. 




Fig. 9 Spatial distributions of the optical/X-ray detected stars in the four age ranges. 



Fig. II lb shows the spatial distribution of all the 8975 optical sources in the Chandra 
ACIS FOV to which no photometric selection was applied. This distribution is clearly 
dominated by field stars and does not show any evidence of a central clustering. The very 
strong absorption of the nebula causes a very patchy spatial distribution of stars. Fig.0 
shows that most of the optical/X-ray sources have age in the range (0.3-10.0) Myr. This 
means that stars with age outside this range are not photometrically cluster members. 
Therefore, we selected all optical stars with age in this range and we defined them as 
"optical candidate cluster members". The spatial distribution of these stars is shown in 
Fig.Ull). Although it is still dominated by field stars, some evidence of star clustering is 
visible. 
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Fig. 10 Age distributions of the optical/X-ray sources in the five distinct spatial subre- 
gions defined in Fig. |5] 

From this sample, we defined as "optical/X-ray cluster members" all the 737 optical 
sources with an X-ray counterpart and age consistent with that of the cluster. The spatial 
distribution of these stars is shown in Fig. II lb . where a strong evidence of clustering is 
seen. The clustered distribution of these stars allows us to confirm the conclusion that 
most of the optical/X-ray detected stars are cluster members, and clearly shows how 
X-ray observations are efficient to select cluster members. We note that these latter stars 
are about 90% of the whole sample of optica l/X-ray sources (828 stars), in agreement 



Damiani et al 



(2004) 



with the field star contamination estimated in 

Fig. Illl i shows the spatial distribution of the optical candidate cluster sample from 
which the 737 optical/X-ray cluster members were subtracted. These stars do not show 
any clustering evidence; on the contrary they show the same patchy pattern of the con- 
taminated sample and this suggests that most of them are field stars. We conclude that 
our optical/X-ray sample suffers from a small contamination, and includes a very large 
fraction of the cluster members within the studied magnitude limit. 



7. Luminosity and Initial Mass Functions 

Fig. 1121 shows the V magnitude distribution of the optical/X-ray cluster members (solid 
histogram) and that of all optical candidate cluster members (dashed histogram) defined 
in Sect. El These distributions can be considered only as lower and upper limits to the 
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Fig. 11 (a): spatial distribution of all the 8975 optical sources within the Chandra ACIS 
FOV; (6): same distribution for the 3427 stars with magnitudes and colors consistent 
with that of the cluster, whose age is between 0.3 and 10.0 Myr (referred to as "optical 
candidate members"); (c): spatial distribution of the 737 Optical/X-ray cluster members; 
(d): spatial distribution of the 2690 optical candidate member sample after subtracting 
the 737 optical/X-ray cluster members. 

Luminosity Function of the cluster. In fact, due to the limited sensitivity of the Chandra 
ACIS detector, we expect that the sample of "optical/X-ray cluster members" is not 
complete but, how we have discussed above, the expected incompleteness should be small 
notwithstanding the sample includes only about 21% of the "optical candidate cluster 
members" . Due to the strong gradient of the spatial density distribution of the whole 
region, it is impossible to choose a suited region of the FOV to be used as "control field" 
to estimate the field star contamination of the optical candidate cluster members. 

In addition, in the case of pre-main sequence stars, the observed luminosity function 
cannot be directly converted into mass function. For these reasons, we have constructed 
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Fig. 12 V magnitude distribution of the X-ray cluster members (solid histogram) and of 
all optical candidate cluster members (dashed histogram) . 



the mass function using the "incomplete" sample of "optical/X-ray cluster members", 
that includes stars with age derived from the CMD between 0.3 and 10 Myr. For these 
stars, we directly estimated the stellar masses from their position on the V vs. V — I 
CMD, as described in Sect. El In order to obtain the mass function of the whole cluster, we 
first considered the completeness of the optical catalog, using the results of the artificial 
star tests, described in Section l3~Tl To obtain the optical complete mass distribution, the 
fractions of the retrieved artificial stars in each V and I magnitude bin, were interpolated 
to the position of the stars in the V vs. V-I CMD. The minimum fraction computed in 
each mass bin is reported in col. 2 of Table El 

To take into account the X-ray incompleteness of member selection, we corrected 
the mass function assuming that the fraction of NGC 6530 members detected as X-ray 
sources is, at the same sensitivity, equal to that detected for a cluster of similar age for 
stars of a given mass. 

The most suitable region for this comparison is certainly the Orion Nebula Cluster, 
one of the best studied clusters, for which all cluster members are known from ex- 
tensive optical and inf rared surveys an d for w hich Chandra X-ray observ ations were 
recently publish e d by 



Flaccomio et al 



Flaccomio et al 



(2003) and 



Feieelson et al 



(2002). Fig. 3 of 



shows the cumulative X-ray Luminosity Functions relative to 



the complete optical sample of stars in the Orion Nebula Cluster, for eight ranges of 
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Table 6. Mass Function of NGC 6530. 



AM Min. optical X-ray fraction for ( aLo^m W" ( ALogM ) ( aLoem ) 1 "" 

(Mq) completeness Log(L x )=29.7 Log(L x )=29.2 
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mass from 0.1 to 50 Mq. In order to obtain an appropriate correction to the NGC 6530 
Mass Function derived from the optical/X-ray star sample, we considered that the sen- 
sitivity of the NGC 6530 Chandra ACIS observation allows us to measure X-ray lumi- 
nosities greater than Log(Lx) — 29.2(erg/sec) in the central region and greater than 
Log(Lx) — 29.7(erg/sec) in the external region of the ACIS FOV for an X-ray spectrum 
with fcT=lkeV, and Nh corresponding to t he optical extinction (Da miani et al., in prepa- 



ration). For each mass range considered by 



Flaccomio et al 



1 2003]), we took the fraction 



of X-ray detected stars above Log(Lx) — 29.2 (erg/sec) and Log(Ljf) = 29.7 (erg/sec), 
given in columns 3 and 4 of Table and we used these fractions to obtain the minimum 
and maximum correction to our incomplete Mass Function. The mass function of 
NGC 6530 AN/ AlogM is therefore taken as the average of the minimum and maximum 
distributions (AN/ AlogM) m i n and (AN/ AlogM) m:ix , resulting from these corrections. 

The results are given in Table El and shown in Fig. 1131 where the dashed histogram 
shows the X-ray incomplete mass function derived from the X-ray cluster members, while 
the black bullets indicate the corrected mass function of the cluster; the error bars show 
the minimum and maximum corrected distributions. Note that, with the exception of 
the lowest mass bin, the optical incompleteness is always smaller than 60% and it makes 
sense to correct it as described before. 

Since the median age of the cluster is smaller than the typical time scale for dynamical 
evolution (~ 100 Myr) and/or stellar evolution, the measured mass function is equal to 
the Initial Mass Function. The cluster IMF increases going down from 6.5 to about 
0.4 Mq, where it shows a peak, and decreases for lower masses. In the mass range (0.6- 



Prisinzano et al.: NGC6530 

x= 1.22±0.17 



29 




1.5 - - 

i ; 
i ; 

l.o L i i I i i i i I i i i i I i ii i 

0.5 0.0 -0.5 

Log(M/M s ) 

Fig. 13 The dashed histogram shows the X-ray incomplete mass function derived from the 
X-ray cluster members, while the black bullets indicate the mass function of the cluster 
corrected according to the procedure described in Sect. (average values between the 
minimum and maximum corrections, indicated in figure by the error bars). The field star 
mass function given in Chabrier (2003), arbitrarily normalized to the present data, is 
plotted as the dotted- dashed line. The solid line is the power law fitting the NGC6530 
mass function in the range (0.6-4.0) Mq. The power law index and the corresponding 
rms uncertainty are also indicated. The vertical dotted line corresponds to the value at 
which the photometric completeness limit is larger than 40%. 



4.0) Mq, it was fitted with a power law of index x — 1.22 ± 0.17 that is consistent with 
th e Salpeter index 1.35. T he lower limit of the mass range was chosen equal to that used 



by 



Muench ct al 



1 2002b) in Orion Nebula cluster, in order to compare the slope of the 



IMF of these similar clusters in the same mass range. The power law index obtained for 



and IC 348 clusters i n the same mass range l|Muench et al 



Muench et al 



Muench et al. 


2002alb 


Briceno et al. 


2002 



2003). In the common mass range, while the IMF of most star forming 



regio ns and open clu sters, such as Trapezi um, IC 348, Lambda Orionis <|Navascues et al 



2QQJJ), and Pleiades <|Bouvier et al 



1998) flattens at about 0.8 Mq, the NGC6530 IMF 



appears to decrease for mass as lower than about 0.4 M©. A quantitative comparison 
with the IMF of other star forming regions or open clusters for very low masses cannot 
be performed with available data, since the NGC6530 IMF is limited to M > 0.2 M . 
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Th e NGC 6530 IMF has also been compared to the field star IMF derived by 



Chabrier 



1 2003(1 that is displayed in Fig. U3 by the dotted- dashed line. The NGC 6530 IMF is 
reminiscent of the lognormal shape derived for the galactic disk field stars. 

The total mass of the optical/X-ray members is about 560 M©; using the corrections 
for optical and X-ray incompleteness we find that the total mass of the cluster, including 
stars down to 0.4 M@, is between 700 and 930 M Q . 



8. Summary and Conclusions 



We used BVI images taken with the WFI camera of the ESO/2.2m, available at the 
ESO/ST-ECF Science Archive, to obtain multi-band photometry of the very young open 
cluster NGC 6530. This cluster is located in a star formation region known as Lagoon 
Nebula (M8) that is one of the brightest nebulae in the solar vicinity. The whole field 
shows evidences of strong obscuration due to the dense molecular cloud located just 
behind the cluster. 

The present photometric catalog reaches down to V ~ 23 and allows us to significantly 
increase the knowledge o f this cluste r . In fa ct, the most recent photometric survey of this 



cluster was published by 



survey 



Sung et al 



Sung et al 



I 2000) and reaches only down to V = 17. From their 
1 2000) derived a cluster distance of 1800 pc corresponding to a distance 
modulus of (V — My) = 11.25, which is much larger than the value d=1250pc, inferred 
from the present work by considering a more complete CMD and taking advantage of 
the total absence of background stars beyond the main sequence locus at the cluster 
distance. 

The near-infrared CMDs of the same FOV, obtained using JHK magnitudes taken 
from the 2MASS catalog, suggest that this region includes a large number of very red- 
dened objects that are not present in the optical CMDs, at least within our limiting 
magnitude. Few of them have been detected as X-ray sources. 

Optical and near-infrared data alone do not allow us to identify the cluster stars in the 
CMDs, since they do not form a well defined sequence but, on the contrary, they populate 
a wide region in the CMDs. For this reason it has been crucial to cross-correlate optical 
data with X -ray sources detected with Chandra ACIS instrument and very recently 



published by 



Da miani et, al 



I 2004). A total of 828 optical/X-ray sources were found and 
placed in the CMD. 90% of these stars are very probable cluster members and identify a 
very well defined pre-main sequence region, that is the cluster locus comprised between 
the 0.3 and 10 Myr isochrones. 

Masses an d age of these clust er members were estimated using the evolutionary tracks 
computed by Siess et al. |200cf) . The age distribution of these stars indicates that the 
median age of the cluster is about 2.3 Myr . A statistically significant trend in the spa- 
tial distribution of these stars was found as a function of age, as already suggested by 
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Lada et al 



(1976J) and 



Damiani et al 



I 2004). This suggests a sequential star formation 



process from north to south. 

The Initial Mass Function derived from the photometrically selected stars with an X- 
ray counterpart was corrected for incompleteness of X-ray data assuming that the fraction 
of NGC 6530 members detected as X-ray sources is, at each given sensitivity, the same of 
that detected in the Orion Nebula Cluster. In the mass range (0.6-4.0) Mq, the corrected 
Initial Mass Function can be represented by a power law with index x = 1.22 ± 0.17, 
consistent with the Salpeter index 1.35, while at smaller masses it shows a peak and then 
it starts to decrease. The resulting IMF is similar to that obtained for coeval clusters, 
such as Trapezium, Taurus and IC348, in the same mass range. 
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